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Abstract 

In this note, we present an eighth-order derivative-free family of iterative methods for nonlinear equations. The pro- 
posed family shows optimal eight-order of convergence in the sense of the Kung and Traub conjecture [5] and is based 
on the Steffensen derivative approximation used in the Newton-method. As a final step, having in mind computational 
purposes, a derivative-free polynomial base interpolation is used in order to get optimal order of convergence with 
only four functional evaluations. Numerical esperiments and few issues are discussed at the end of this note. 
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1. Introduction 

Let / : D c % — > 21 be a sufficiently differentiable function of single variable in some neighborhood D of a, 
where or is a simple root (f'(a) + 0) of nonlinear algebraic equation f(x) = 0. The well-known Newton method is 
denined by the iteration 



x n +\ 



f(Xn) 



(1) 



/'(*„) 

which shows a second-order convergence. One can easily get Steffensen's approximation for first order derivative as 



f(x„ - Kf(x„J) ~ f(x„) - Kf(x„)f'(x„), 

Kf(x„)f'(x n ) * f(x n ) - f(x„ - Kf(x n )), 

^ 1 /(■*«) ~ f(Xn ~ Kf(x„)) 
K f(Xn) 



(2) 



f'(x n ) 



If we substitute the derivative approximation © in ([T}, we obtain Steffensen's second order accurate derivative-free 
iterative method for non-linear equations yl . 



— X n Kf{x n ), 

f(x„) 2 



X n +\ — X n — K 



(3) 



f(Xn) ~ f{Wn) 
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In 2012, an optimal eighth-order iterative method [2] was proposed by Y. Khan et al. as 

f(x„) 



Zn 



f(x n y 

d fiyn) ] f(y n ) 

fJ- f(Zn) 



fi + vq 2 „f'(z„y 

where /j. + 0, v € 21, q n = f(z„)/f(x n ), G(t) is a real-valued function with G(0) = 1, G'(0) = 2, G"(0) < oo, and 

f'(z„) ~K-C{y„-Zn)-D{y n -z n ) 2 , 
f(x n ) - f(y n ) 



/fa) - fjZn) 

y n -z n ' ( 5 ) 
f\x n )-H H-K 



(%n yri)(Xn Zn) (?Cn Zn)^ 



H 
K 
D 

{%n ~ yn 

H-K 

C = — D(x n +y n -2z„). 

Xn Zn 

In the original draft of paper [2] the expression for C has typo-mistake, which is corrected here. Actually (|5) poly- 
nomial interpolation approximation for f'(z n ) is given in y|]. Clearly (|4]i iterative scheme is not derivative free. The 
main contribution in this paper is to use the idea of iterative scheme <j4j by introducing Steffensen's derivative ap- 
proximation for f'(x n ) and then finally construct derivative-free approximation for f'(z n ) without reducing order of 
convergence. 



2. Construction of derivative-free family 

First we construct an interpolation polynomial approximation for /'(z„). Suppose we have f{x n ), f(w„) (defined 
in[3]l, f(y„) and f(z n ), One could construct a three-degree polynomial as follows 

(p(</>) = f(y„) + ri(0 - y„) + r 2 (0 - y„f + r 3 (cf> - y n f, 
\p'(<f>) =n +2r 2 (0-y„) + 3r 3 (0-y„) 2 . 

By using four functional values, we get the following system of equations: 

2 3 

v\ — r\a + r 2 a + r 3 a , 

v 2 = rib + r 2 b 2 + r 3 b 3 , (7) 

2 3 

V3 = r\c + r 2 c + r-$c , 
where 



(8) 



After solving (0 for r\, r 2 and r 3 and substituting them in (O implies the following approximation for/'(z„): 

, b(b - c) vi -3b 2 + 2bc + 2ab - ac v 2 b(b - a) v 3 

/ (z„) « ij/ n = — -— + — — + - — -— . (9) 

(a — b){a — c) a (a — b)(b — c) b (a - c)(b - c) c 

2 



Vl 


= f(x„) - 


f(jn), 


V'2 


= f(z„) - 


f(y n ), 


V 3 


= f(w„) - 


-f(y n ), 


a 


— x n — y n 




b 


= z„-y„ 




c 


= w„ - y„ 





We consider the following family of iterative methods for nonlinear equations: 

w„ = X„- Kf(x n ), 

f(Xn) 2 



f(x„)-f(w n ) 

f(yn)fUn) r(f (10) 

z» -y n - k— — - — - — -G(h,t 2 ), 

f(x„) - f(w„) 

X n +\ = Z n ; H(SUS2), 

. , f(y n ) t f(y„) f(z n ) f(z„) „ 

where h = -7—, h = — — -, si = 7-—, s 2 = — — - and *<gfc 0) e !R. 
/(*„) /(w„) /(*„) /(w„) 

3. Convergence analysis 

We state the following theorem about the order of convergence of the family described in (TTOb . 

Theorem 1. Let f : D c % — > S X foe a sufficiently differentiable function, and a e D is a simple root of f(x) — 0, for 
an open interval D. If xq is chosen sufficiently close to a, then the iterative scheme given in ( 1701 ) converges to a. If G 
and H satisfy 



G(0,0)= 1, — 

Ot\ 



= 1-? 
dh 

(0,0) 2 



= 1,H(0,0)=1, ^ 

(0,0) 



ds\ 



(0,0) 



-*ir 

OS2 



= 0, (11) 

(0,0) 



d^G d^G d^G cP'H d^H d^H 

and — — , — — , , — — , — — , are bounded at (0,0) then the iterative scheme ( 17 01 ) shows an order of 

dt\ dt\ dhdt 2 ds\ dsj dsids 2 [] —" 

convergence at least equal to eight. 



1 f (k \a) 

Proof. Let the error at step n be denoted by e n — x n — a and let us define c\ - f'(a) and c* = , k = 2,3, • • •. 

k\ f'(a) 

If we expand / around the root a and express it in terms of powers of error e„, we obtain 

f{x n ) = c\(e n + c 2 e\ + c^e\ + c A e A n + c 5 e^, + c b e\ + c-,e n + c%e\ + 0{e 9 n )), (12) 

9 9 9 9 999 99 

/(w„) = — Ci(— 1 + KC\)e n + c\c 2 (—2>kc\ + I + k c l )e n — c\{Akc\Ct, + 2c 2 kc\ — 2k c l c 2 — C3 - 3c3^C[ 

-IT'} 99 9 9 ^ ^ 99 ^ 

+ 6c 4 jTCi - 4c 4 /c J Cj (13) 

+ C4* 4 cM + -.- + 0(«2), 

9 999 9999^ 

y„ — a = —C2(— 1 + Kc{)e n + (2c3 — 3a"CiC3 + C3K Cj + 2c 2 /fci - 2c 2 - a: c 1 c 2 )e n + (3c 4 + 10c2AriC3 — 6kc\C/[ 

, A 2 2 3 3 t 2 2 -7 c 3 , o 33 , o 2 2 3 , ., 3 3 3 3\ 4 

+ 4-C4K C[ - C4/f Cj - 7 a: C[C2C3 - /QC3 - 5atic 2 + 2c2^ C[C3 + 3a: CjC 2 + 4c 2 - c l c 2 )e n 
+ ■■■ + 0(e 9 n ), 

(14) 

9 999 9999"^ 

- — ciC2(— 1 + KC\)e n - c\{— 2C3 + 3/cciC3 - C3/c C[ - 2c 2 /fCi + 2c 2 + a- C[C 2 )e K - ci(-3c4 - 10c2a:ciC3 

99 "X ~\ 99 — T 99"^ "X 

+ 6kc\Ch — Ac/[K Cj + C4/C Cj + Ik c x C2C^ + 7c2C3 + 1 kc\c 2 — 2c 2 icc\c^ — 4k c x c 2 — 5c 2 
+ K i c\cl)e A n + --- + 0{el), 



f(y n ) 

f(Xn) 

f(y«) 

f(w„) 



(15) 

99 9 99999 Q 

= —c 2 (— 1 + KC\)e n + (2c3 - 3a:ciC3 + C3/C Cj + 3c 2 a-ci - 3c 2 - k c l c 2 )e n + ■ • ■ + 0(e n ), (16) 
: c 2 e„ + (-kc\Cj + 2c\kc\ + 2c3 - 3c 2 )ef 7 + ■ • • + 0(el). (17) 



3 



The Taylor series expansion of G(t\ , £2) is given by 

2 2 

J f(y n ) f(y n ) \ , , /(y«) , /(y«) , , / /(?«) \ , / /(yJ \ , , ( f(y n ) \( f(y„) \ , n / 3 j\ n ^ 

By using O, CQh CGK GSl, O, we find 

Z n — os = (— A\c\ + 6k 2 c\c\ — C3 + 2a:ciC3 - IOcjKCi - C3/c 2 C[ + 5c 2 + Aic 2 kc\ + 2c\A 3 kc\ — c 2 2 A 3 k 2 c\ 
+ 3A 2 c\kc\ - 3A 2 c\k 2 c\ + A 2 c\ic'c\ - A 2 c\ - k 3 c\c 2 2 - c^A-s)c 2 e A n + {-Ai^c\c\ - 4k 2 c\c4C 2 

- 3\A 2 c 4 2 kc x + 36A 2 c 4 2 k 2 c\ - \9A 2 c\i^c\ + 4A 2 c 4 2 k 4 c\ - 23c 4 2 A 3 kc x + ISc^k 2 ^ - 5c 4 2 A^c\ 

1 422 442 3 3 222 332 

- \5A\c 2 kc\ + 6A\c 2 k Cj + 5c4KCic 2 + 3k c x c 2 c 3 + ca,k c x c 2 + 68k c x c 2 c 3 — 25c 3 k ~ c { c 2 — 2c^c 2 

- 1%kc\c\c 3 — 3k 4 c\c 4 2 + 5kc\c\ + c\k 3 c\ — 2c\ — 36c 4 + 32c\c 3 — 66c\k 2 c\ + %0c 2 kc\ + 24k 3 c\c 4 2 
+ 9A\c\kc\c 3 - 3A\c\cj,k 2 c\ + 2\A 2 c\kc\c 3 - 21A 2 c\c 3 k 2 c\ + \5A 2 c\c 3 K i c\ - 3A 2 cIk 4 c 4 C3 

+ \5c\A-skc\c?, - \2c 2 2 A 3 c 3 k 2 c\ + 3c 2 2 A 3 c 3 k 3 c\ - 6A\c 2 2 c 3 - 6A 2 c\c 3 - 6c\A 3 c 3 + lOAiCj + IOC2A3 
+ \QA 2 c 4 2 )el + --- + 0{el), (19) 





f{Zn) = C\{—A\c 2 + 6k c x c 2 - C3 + 2kc\Cj, - \Qc 2 kc\ - cjk c x + 5c 2 + A\c 2 kc\ + 2c 2 A 3 kc\ — c 2 A 3 k c x 
+ 3A 2 c 2 2 kc\ - 3A 2 c\k 2 c\ + A 2 c\k^c\ - A 2 c\ - k 3 c\c\ - c 2 2 Ai)c 2 e 4 n + c\{-4i^c\c\ - 4k 2 c\c4C 2 

- 3\A 2 c 4 2 kc x + 36A 2 c 4 2 k 2 c\ - \9A 2 c 4 2 k 3 c\ + 4A 2 c 4 2 k 4 c\ - 23c 4 2 A 3 kc x + \%c 4 2 A 3 k 2 c\ - 5c 4 2 A 3 k 3 c] 

1 422 442 3 3 222 3 32 

- \5A\c 2 kc\ + 6A\c 2 k Cj + 5c/[KC\c 2 + 3k c y c 2 c 3 + C4K c x c 2 + 68k c y c 2 c 3 — 25c 3 ic c\c 2 — 2c\c 2 

- 1%kc iC 2 2 c 3 - 3k 4 c\c\ + 5kc iC \ + c 2 3 k 3 c\ - 2c\ - 36c\ + 32c 2 2 c 3 - 66c 4 2 k 2 c\ + 80c^ Cl + 2\k 3 c\c\ ^ 

2 222 2 222 23 3 24' 

+ 9A\c 2 kc\c 3 — 3A\c 2 c 3 k Cj + 2\A 2 c 2 kc\Ct, — 21A 2 c 2 c 3 k c { + \5A 2 c 2 c 3 k c\ — 3A 2 c 2 k 

% 'X 

+ \5c 2 A 3 KC\c-i — \2c 2 A 3 c 3 k Cj + 3c 2 A 3 c 3 k c\ — 6A\c 2 c 3 — 6A 2 c 2 c 3 — 6c 2 A 3 c 3 
+ 10A lC ^ + 10c^A 3 + lOA 2 c 4 2 )el + ■■■ + 0(e 9 n ), 
ip„ = (2A 2 c\k 3 c\ - 2k 3 c\c\ + 12k 2 c 2 c 3 2 - 2k 2 c\c 2 c 3 - 6A 2 c\k 2 c\ + c 4 k 2 c 2 - 2c 3 2 A 3 k 2 c\ - 20kc\c\ 

+ 6A 2 c 2 kc\ — 2kc\Ch + 2A\c 2 kc\ + Ac 2 kc\c 3 + Ac 2 A 3 kc\ + I0c 2 — 2A\c 2 - 2c 2 A 3 - 2c 2 c 3 (21) 

- 2A 2 c^ + c 4 )c x c 2 e 4 + ■■■ + 0(e 9 n ). 

Finally, H(s\, s 2 ) has Taylor's expansion 

\f(Xn) f(Wn)J \f(Xn)J \f(w n )) \f(-X n )J\f(w n )J KV2 ' 

From (120b and (l2TT i. we deduced the following error equation which leads to the desired result 

e n+l = c 2 2 (-6k 2 c\c 3 c\ - IOAiCj - \QA 2 c\ - IOC2A3 + 6c 2 k 2 c\c\ + 3\k 2 c\c 4 c\ - Ac^c 3 ^ 

33 3 33 3 3 2 !42 I 1 2 2 

+ 4c 3 kciC4 + 4c 3 C4K c\ + \6c 3 ic c x c 2 — 23c\k c\c 2 + 8C4K c x c 2 — C4K — 4c 2 kc\c 3 — 2Qkc\C4C 2 

- C3C4 - 62c 3 k 2 c\c\ + k 4 c\c 2 c\ — 16k 4 c 4 c 2 ci - \2k s c\c\ + 2k 5 c\c 3 c\ - k 5 c\c4c\ + 6A\c 3 2 k 2 c\c 3 

'X "\ "X "X "X "X "X "i "X 

+ 3A\c 2 kc\C4 - 6A\c 2 kc\Cj, - 2A\C 2 k c\c 3 — 3A\c 2 C4icc\ + A\c 2 C4K c\ - 8c 2 A 3 k c\cj, + 4c 2 A 3 kc\C4 

"X "X "X "X 

— 6c 2 A 3 k c x C4 + \2c 2 A 3 k C[C3 + 4Qc 2 c 3 kc\ — A\c 2 C4 + 2A\c 2 c 3 — A 2 c 2 C4 +2A 2 c 2 c 3 — c 2 A 3 C4 
+ 2c 2 A 3 c 3 + 25c 2 + 2c 2 A 3 k 4 c 4 c 3 + 4c 2 A 3 k 3 c 3 C4 — c\A 3 k 4 c\c4 — %c 3 2 A 3 kc\c 3 — 20A 2 c 2 k 3 c 3 c 3 

+ 5A 2 c\kc\C4 - \0A 2 c\k 2 c\c4 + 20 A 2 c\k 2 c\c 3 + \0A 2 c 3 2 k a 'CjC 3 + \0A 2 c\k 3 c\c4 - 5A 2 c\k 4 c\c4 

4 



- \0A 2 c\kc\Ct, - 2A 2 cl^c 5 l C3 + A 2 c\k 5 c\c4 + c 2 c\ - \Oc\ci + 5c\ca + \60k 2 c\c\ - \?>Oic'c\c\ 

- \00c\kc\ + 56k a c\c\ - 32Aic\k 2 c\ + 30Aic 5 2 kci + \AAic\ic'c\ + A6c\A- i ic'c\ - 62c 5 2 A 3 k 2 c 2 

- 16c 5 2 A 3 k 4 c 4 + 40c 5 2 A 3 kci + 108A 2 c^ 3 <^ - 102A 2 c 5 2 k 2 c 2 - 62A 2 c\k 4 c\ + 5QA 2 c 5 2 kci + 18A 2 c^ 5 c^ 
+ A\c\ + A 2 2 c\ + c 5 2 A\ - 2A\c\kc\ - 2A\c\k 4 c\ + A\c s 2 k 2 c\ - Ac\A\k*c\ + dc\A\i?c\ - Ac\A\kc\ 

+ 2c 5 2 A i K 5 c\ + c 5 2 A\k 4 c\ + 15A 2 2 c 5 2 k 4 c 4 - 20 A 2 2 c 5 2 k 3 c\ + \5A 2 2 c\k 2 c\ - 6A 2 2 c\k Ci - 6A\c? 2 ^c\ 

+ A 2 2 c\i£c\ - 2A 2 c\i£c\ - %A\c\A 2 K i c\ + \2A\c 5 2 A 2 k 2 c\ + 6A\c 5 2 A 3 k 2 c 2 - %A\c\A 2 kc\ - 6A\c 5 2 A 3 kc\ 

+ 2A x c\k a c\A 2 - 2A 1 c 5 2 k 3 c]Ai + \Qc\A 3 k a c\A 2 - 20^A 3 ^c 3 A 2 + 20c 5 2 A 3 k 2 c 2 A 2 - 10c 5 2 A 3 k Ci A 2 

- 2c 5 2 A 3 k 5 c 5 1 A 2 + 2A lC ^A 3 + 2A lC ^A 2 + 2A 2 c\A 3 + K 6 c\c 5 2 )el + 0(e 9 n ). (23) 

□ 

It is clear that the considered family of numerical schemes requires four functional evaluations and attains optimal 
convergence order eight according to Kung and Traub conjecture which can be stated as follows [5]: if n is the 
total number of functional evaluations per iteration, then the optimal convergence order of the associated numerical 
procedure is 2 . 



4. Numerical Results 

r — l 

Definition 1. The computational order of convergence Iffl, can be approximated by 

coc „ fa|Q«+i -a){x n -aY l \ ^ 
ln\(x n -a)(x„_i -or) -1 1' 
where X n -\, x„ and x„ + i are successive iterations closer to the root a of f(x) = 0. 

For the purpose of comparison between newly developed family and other derivative-free methods, a list of 
derivative-free methods for nonlinear equations is presented here. 

4.1. The Kung-Traub Eighth-order Derivative-free Method (K-T) 

The Kung-Traub eighth-order derivative-free method is discussed in f^Q], and also considered in [7] is given as 



w„ 

y n 

Zn 

x n +i 



= x n +Pf(x„) 
= x, 



= y n - 







f(Wn)-f(x„)j 




f(x n )f(w n ) \ 


1 


f(yn)-f(x„)j 


f[w n ,x n ] 


f(w„)f(x n )f(y„)\ 



(25) 



f(z„) - f(x„) 
1 



f{Zn)~ f{w„) 



1 



1 



f[y n ,z„] f[w„,y n ] 



1 



f(y n ) - f{x„) 



l 



l 



f[w„,y n ] f[w„,x n ] 



4.2. R. Thukral Ml, M2, M3 Methods 

In 2011, R. Thukral [7] presented three variants of his proposed eighth-order three-point derivative-free method. 
Three members of the family called by author namely, Ml, M2, and M3, are listed as 

^ / ( v.. I i ' 



01 

02 



f(y») 



l + 



f(w„) 



f(y n ) , / f(y„) \2 



5 3 = 



f{Wn) VOn)' 

f[x„,w n ] 

f[w n ,y n Y 



(26) 

(27) 
(28) 



and 



w„ 



x n+l 



= x n +/3f(x n ), 

Pf(Xn? 



y n -4>k 

f. 



/W - /(*«) 



f[Xn,y n ] 
f(Zn) 



(29) 



f(w„) 2 f(x n ) )\f[y„, z„]f[x n , z n ] 



where = 1,2,3, /? E 2V*", 0* are listed in (26)-(28). (|29]l is called Ml, Ml and M3 for 0i, 2 and respectively. 

4.3. Petkovic et al. Type Methods 

In J3l, author developed Petkovic type 1 (PI) and type 2 (P2) derivative-free methods for the comparison of 
numerical efficiency, (PI) and (P2) respectively, are written as 

w n = x„ +/}f{x n ), 

Pf(Xn) 2 



z„ = y n - 1 1 + 



x n+ 1 



f(Wn)-f(x„))' 

f(y n ) , /Oh) 



f(w„) f(x„) 

f{ZnV' ] 



(yv n - x n )f(y„) 



f(Wn) - f(Xn) 



and 



1 - 



w n = x„ + 
y n — x n 



fijnf 



fiyn) \ \l f[x„,y n ]f(Zn) 



2f(y n ) 3 

f(w„) 2 f(x n ) f(w n )f(x„) 2 \f(w n )j j\f[y„,z n ]f[x„,z„]J' 



Pf(x n ), 



■y n 



f(Wn)-f(x„)J' 

l + f(y n )f(x„y l \( f(y n )(w n - X„) 



i-f(yn)f(w n y l !\f(w n )-f(x n ) 



Xn+l 



f(Zn) 

f(w„) 



2f(y n y 



f(y n y 



f[x„,y„]f(z„) 

f(Wn) 2 f(x„) f(Wn)f(x n ) 2 )\f\y„, Z n ]f[x„, Z„] 



4.4. Proposed family (L) 

We define the following weight functions: 



Gi(fi,f 2 ) = 



1 



i - (?i + r 2 ) + «x?i + t 2 y 



G 2 (fi,f 2 ) = l + h+t 2 + % + 1.9f 2 +4.4fjf 2 , 

ffl(*l,J>z) = 1, 

1 

Hl{s ] _,S 2 ) 



1 + S\S2 + s 2 + s\ ' 



H 3 (s u s 2 ) = 1 + s 2 + s 2 , 
Hi(si, S2) - I + s\ + si + 2s\S2, 



H 5 (si,s 2 ) = 



1 

1 — 20sis 2 ' 



(30) 



(31) 



(32) 

(33) 
(34) 

(35) 

(36) 
(37) 

(38) 



Functions 



Roots 



fi(x) = exp(x) sin(x) + ln{\ + x 2 ) a = 

f 2 (x) = x 15 + x 4 +4x 2 - 15 a = 1.148538... 

f 3 (x) = (x- 2)(x 10 + x + 1) exp(-x - 1) a = 2 

fa{x) = exp(-x 2 + x + 2) - cos(x + 1) + x 3 + 1 a — -1 

/sW = ( x + 1) exp{sin{x)) - x 2 exp{cos(x)) - 1 a = 

/ 6 (x) = sifi(i) 2 - x 2 + 1 a = 1.40449165... 

/ 7 (x) = 10 exp(-x 2 ) - 1 or = 1.517427... 

Mx) = (x 2 - - 1 a = 1.414214... 

/ 9 (x) = /«(x 2 + x + 2)-x+l a = 4.15259074... 

f w (x) = cos(x) 2 - x/5 a = 1.08598268... 

x 

f n (x) = sin(x) - - a = 

f u (x) = x 10 - 2x 3 - x + 1 a = 0.591448093... 

/i 3 (x) = exp(sin(x)) - x + 1 a = 2.63066415... 



(/n(*),*0 ) 


L 


K-T 


Ml 


M2 


M3 


PI 


P2 


/i,0.25 


(LI) 6.38e-247 


3.14e-136 


1.69e-141 


7.43e-142 


1.69e-141 


3.20e-113 


8.98e-120 


/2, LI 


(LI) 1.2376e-652 


3.72e-61 


3.44e-62 


3.44e-62 


3.44e-62 


2.68e-7 


2.67e-7 


/3.2.1 


(LI) 1.057e-422 


1.91e-60 


1.49e-60 


1.49e-60 


1.49e-60 


7.71e-8 


7.56e-8 


/ 4 , -0.5 


(LI) 2.952e-383 


5.11e-362 


1.92e-362 


1.93e-362 


1.92e-362 


9.99e-367 


8.78e-366 


/s, 0.25 


(LI) 2.336e-407 


4.13e-328 


6.52e-326 


9.47e-326 


6.52e-326 


1.98e-322 


2.56e-332 


/e, 1.2 


(L8) 1.719e-421 


1.00e-327 


4.58e-341 


7.57e-344 


4.58e-341 


1.79e-381 


1.72e-405 


/7,2 


(L2) 7.264e-238 


5.19e-88 


1.24e-120 


6.40e-124 


1.24e-120 


1.51e-187 


6.79e-228 


/s, 1.7 


(L3) 1.429e-234 


1.23e-113 


1.74e-171 


5.45e-188 


1.74e-171 


5.96e-211 


4.84e-167 


/9,4.4 


(L4) 2.504e-997 


1.15e-928 


4.52e-942 


1.27e-965 


4.52e-941 


6.15e-904 


4.11e-937 


/io, 1-5 


(L5) 2.81e-305 


7.19e-303 


5.07e-284 


1.84e-245 


5.07e-285 


4.91e244 


1.78e-275 


/i i,0.25 


(L6) 2.35e-1143 


3.65e-782 


1.00e-819 


4.98e-823 


1.00e-819 


5.13e-794 


5.13e-812 


/i2, 0.25 


(L6)7.86e-318 


2.03e-256 


5.65e256 


1.82e-254 


5.65e-256 


1.07e-264 


6.31e-268 


/is, 2.0 


(L7) 2.54e-436 


2.63e-396 


1.94e-378 


5.1e-378 


1.94e-378 


8.70e-380 


6.80e-379 


(COC) 


/i 


(LI) 7.9999 


7.9986 


7.9995 


7.9998 


7.9995 


7.9958 


7.9978 


h 


(LI) 8.0000 


7.8671 


7.9371 


7.9371 


7.9371 


3.2715 


3.2731 


h 


(LI) 7.9999 


7.8660 


7.9047 


7.9047 


7.9047 


4.2595 


4.2675 


U 


(LI) 8.0000 


7.9905 


7.9905 


7.9905 


7.9905 


7.9907 


7.9907 


h 


(LI) 8.0000 


7.9884 


7.9882 


7.9882 


7.9882 


8.0000 


8.0000 


k 


(L8) 8.0000 


8.000 


8.0000 


8.0000 


8.0000 


8.0000 


8.0000 


fi 


(L2) 7.9999 


8.0097 


8.0025 


8.0018 


8.0025 


8.0005 


8.0001 


h 


(L3) 8.0000 


8.0027 


8.0004 


8.0002 


8.0004 


8.0000 


8.0002 


h 


(L4) 8.0000 


8.0000 


8.0000 


8.0000 


8.0000 


8.0000 


8.0000 


fw 


(L5) 8.0000 


8.0002 


7.9845 


7.9797 


7.9845 


8.0003 


7.9833 


fn 


(L6) 11.000 


10.996 


10.996 


10.996 


10.996 


10.996 


10.996 


hi 


(L6) 8.0000 


8.0000 


7.9809 


7.9807 


7.9809 


7.9822 


8.0000 


fn 


(L7) 7.9999 


8.0000 


8.0000 


8.0000 


8.0000 


8.0000 


8.0000 



Table 1: Numerical comparison between three-point derivative-free methods 

7 



where f, and s, are defined in ([Tol l. Further we give names to methods for the purpose of simplicity as follows 

(Gi, H u a) = +0.01, k = 0.01), L2 = (G u H u u = -0.022, k = 0.01), 
(d, Hi, w = -0.001, k = 0.01), LA = (G 2 , Hi, w = +0.01, k = 0.01), 

(Gi, H 3 , u) = -0.01, k = 0.01), L6 = (G,, // 2 , w = +0.01, k = 0.01), <J ' 

(Gi, tf 4 , w = +0.01, k = 0.01), L8 = (G,, H 5 , u = +0.01, k = 0.01). 

A set of thirteen nonlinear equations is used for numerical computations from [7], in Table Q] All the families in 
the numerical implementation are derivative-free and use four function evaluations to get the order of convergence 
eight . For all methods, 12 (TNFE) total number of function evaluations are used, and absolute error (\x n - a\) is 
displayed. Computational order of convergence is calculated according to (l24l > for the method. All numerical values 
for methods K-T, Ml, M2, M3, PI, P2 are taken from 0]. 

5. Conclusion 

In this note, we have presented a family of eighth-order derivative-free methods. The proper selection of weight 
functions showed a reasonable reduction in error as compared to other referenced derivative-free methods. It is 
obvious that constructed family has broad choice for the weight function in the third and fourth step of the method. 
The true essence of the family is hidden in the construction of interpolation polynomial for the approximation of f'(z) 
and weight functions make it more flexible to get higher performance and efficiency. 
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